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ABSTRACT 

I present a simplified analytical model that simulates the evolution of the binary pop- 
ulation in a dynamically evolving globular cluster. A number of simulations have been 
run spanning a wide range in initial cluster and environmental conditions by taking 
into account the main mechanisms of formation and destruction of binary systems. 
Following this approach, I investigate the evolution of the fraction, the radial distribu- 
tion, the distribution of mass ratios and periods of the binary population. According 
to these simulations, the fraction of surviving binaries appears to be dominated by 
the processes of binary ionization and evaporation. In particular, the frequency of bi- 
nary systems changes by a factor 1-5 depending on the initial conditions and on the 
assumed initial distribution of periods. The comparison with the existing estimates 
of binary fractions in Galactic globular clusters suggests that significant variations in 
the initial binary content could exist among the analysed globular cluster. This model 
has been also used to explain the observed discrepancy found between the most recent 
N-body and Monte Carlo simulations in the literature. 

Key words: stellar dynamics ~ methods: analytical - stars: binaries: general - 
Galaxy: globular clusters: general 



1 INTRODUCTION 

The study of the evolution of binary stars both in the solar 
neighborood and in stellar clusters is one of the most inter- 
esting topic of stellar astrophysics. Binarity, under particular 
conditions, induces the onset of nuclear reactions leading to 
the formation of peculiar objects like novae and determines 
the fate of low-mass stars leading to SN la explosions. In col- 
lisional systems binaries provide the gravitational fuel that 
can delay and eventually stop and reverse the process of core 
collapse in globular clusters (see Hut et al. 1992 and refer- 
ences therein). Furthermore, the evolution of binaries in star 
clusters can produce peculiar stellar objects of astrophysic 
interest like blue stragglers, cataclysmic variables, low-mass 
X-ray binaries, millisecond pulsars, etc. (see Bailyn 1995 and 
reference therein). 

From an observational point of view, there is an exten- 
sive literature on the analysis of the main characteristics of 
the binary population of the Galactic field (Duquennoy & 
Mayor 1991; Halbwachs et al. 2003 and references therein). 
Conversely, the analysis of the binary population in star 
clusters is still limited to small samples of open clusters 
(Bica & Bonatto 2005) and to few globular clusters (Yan 
& Mateo 1984; Romani & Weinberg 1991; Bolte 1992; Yan 
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& Cohen 1996; Yan & Reid 1996; Rubenstein & Bailyn 1997; 
Albrow et al. 2001; BeUazzini et al. 2002; Clark, Sandquist 
& Bohe 2004; Zhao & Bailyn 2005). More recently, Solhma 
et al. (2007) estimated the binary fraction in 13 low-density 
Galactic globular clusters. They found that the fractions of 
binary systems span a wide range comprised between 10- 
50%. This spread could be due either to differences in the 
primordial binary fractions or to a different evolution of the 
binary populations. 

Theoretical studies addressed to the study of the evo- 
lution of the properties of the binary population in globular 
clusters are mainly based on two different approaches: i) 
full N-body simulations and ii) Monte Carlo simulations. 
The former follows the dynamical evolution of the stellar 
system assuming simplified treatments of single and binary 
star evolution (see Portegies Zwart et al. 2001; Shara & Hur- 
ley 2002; Hurley & Shara 2003; Trenti, Heggie & Hut 2007; 
Hurley et al. 2007). This approach is extremely expensive 
computationally and have been often performed with un- 
realistically small numbers of binaries. The latter method 
uses a binary population synthesis code to evolve large num- 
bers of stars and binaries by introducing a simple treatment 
of dynamics. In this type of approach it is often assumed 
that all the relevant parameters of the cluster (central den- 
sity, velocity dispersion, mass, etc.) remain constant during 
the cluster evolution (see Hut, McMillan & Romani 1992; 
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Di Stefano & Rappaport 1994; Sigurdsson & Phinney 1995; 
Portegies Zwart et al. 1997; Portegies Zwart, Hut & Verbunt 
1997; Davies 1995; Davies & Benz 1995; Davies 1997; Ra- 
sio et al. 2000; Smith & Bonnell 2001; Ivanova et al. 2005). 
The most recent works performed following these two differ- 
ent approaches leaded to apparently contradictory results: 
while Monte Carlo simulations (Ivanova et al. 2005) predict 
a strong depletion of the binary population in the cluster 
core, the N-body simulations performed by Hurley et al. 
2007 show a progressive increase of the core binary fraction 
as a function of time. 

In this paper I study the evolution of the binary pop- 
ulation in a dynamically evolving globular cluster using a 
simplified analytical approach. Following this approach, I 
investigate the evolution of the fraction, the radial distri- 
bution, the distribution of mass ratios and periods of the 
binary population. A number of simulations have been run 
spanning a wide range in initial cluster and environmental 
conditions. 

In Sect. [2] I describe in detail the code. In Sect. |3] I 
present the complete set of simulations and discuss the de- 
pendence of the results on the underlying assumptions. In 
Sect. |4] the predicted evolution of the binary fraction as a 
function of the initial cluster structural and environmental 
parameters is presented. Sect. [5] is devoted to the study of 
the radial distribution of binary systems. In Sect. [6] and [7] 
the distribution of periods and mass ratios are investigated 
together with their dependence on the environmental condi- 
tions. A comparison with the most recent N-body and Monte 
Carlo simulations is performed in Sect.|8l In Sect. |9] the pre- 
dictions of the code are compared with the estimates of bi- 
nary fractions in globular clusters available in the literature. 
Finally, I discuss the obtained results in Sect. 1101 

For clearity, I list below the notation used throughout 
the paper. 



General 




Pr 


density at the distance r from the 




cluster center 




velocity dispersion at the distance 


Tc 


r 

core radius 


rt 


tidal radius 




binary fraction 




core binary fraction 


V 


fraction of evaporating systems 


t9 


cluster age 


a 


IMF power-law index 


<m> 


average object mass 




(binaries-|-single stars) 


sys 


number of cluster objects 




(binaries-|-single stars) 


M 


cluster mass 


Mc 


cluster core mass 




local relaxation time at the 




distance r 




relaxation radius 


X 


polytropic index 



G 

log A 
-I>r 



Binaries 

mi 

m2 

P 

V 

Nb(mi,m2, log P) 
nb.r{mi,m2, log P) 



g{log P) 
1 

Smgle stars 
ma 

Ns{ms) 



ns,r{ms) 



distribution of velocities at the 
distance r 

cross section of the process X 
escape velocity at the distance r 
Newton constant 
Coulomb logarithm 
potential at the distance r 
energy per unity of mass at the 
distance r 

mass of the primary star 
mass of the secondary star 
period 

hardness parameter 
binding energy 
number of binaries with 
parameters mi,m2,log P 
number density of binaries with 
parameters mi,m2,log P at the 
distance r 

distribution of periods 
mass ratio 

mass 

Number of single stars with mass 
nis 

number density of single stars 
with mass nia at the distance r 



2 METHOD 

Each simulation is characterized by four initial parameters: 
the cluster central density po, the central velocity dispersion 
(T„,o, the initial binary fraction ^ and the rate of evaporating 
systems v. 

The code simulates the evolution of the number of bi- 
naries and single stars by taking into account the main pro- 
cesses of formation, destruction and evolution of binaries 
and single stars. In particular, the following processes have 
been considered: 

• tidal capture 

• direct collisions 

• exchanges 

• collisional hardening 

• binaries ionization 

• stellar evolution 

• coalescence 

• mass segregation 

• evaporation 

A detailed description of the analytical treatment of 
each of these processes is given in the sub-sections from 12.21 
to l2.7l Moreover, during its evolution, the cluster dynamical 
parameters are assumed to evolve (see Sect. I2.8"1l . 

The binary population of the system has been divided 
in several groups according to the masses of the primary 
and secondary stars (mi and m2), and to the period (P). 
Similarly, the population of single stars has been divided in 
groups of masses (ms). 

The evolution of each population of single and binary 
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stars is calculated in layers of variable width located at dif- 
ferent distances from the cluster center. 

At the beginning of the simulation, the density and the 
distribution of velocities of both single and binary stars is as- 
sumed to follow the radial behaviour of a mono-mass King 
profile, as expected for a non-relaxed stellar system. Dur- 
ing this initial stage, the cluster is assumed to be a non- 
coUisional system. Therefore, the only processes at work are 
those related to stellar evolution and evaporation. 

After this initial stage, the evolution of the cluster and 
its binary population is followed in time-steps. For each 
time-step the adopted general procedure is schematically the 
following: 

(i) The code calculates the maximum radius in which 
the local relaxation time is smaller than the cluster age 
{relaxation radius; rrei)- The local relaxation time at the 
radius r has been calculated using the relation 



tr = 0.34 



G^pr < rn> log A. 



Binney & Tremaine (1987) (1) 



with A = OANsys- 

(ii) The number of newly formed and destroyed bi- 
nary systems have been calculated in all the layers lo- 
cated inside the relaxation radius. The number of binaries 
Nb{mi,m2,log P) and single stars A'^s('tIs) belonging to each 
group have been consequently updated; 

(iii) A set of partially-relaxed multi-mass King profiles 
has been calculated (see Appendix A). The density and ve- 
locity distribution of each sub-population of binary and sin- 
gle stars have been assumed to follow a given profile accord- 
ing to their systemic masses; 

(iv) The mass of the cluster has been calculated and the 
values of po, cr„,o and Wo have been updated (see Sect. 12.8)) : 

Point i) to iv) have been repeated until the cluster reach 
an age of 13 Gyr. In the next sub-section I will report the 
assumptions made in all the simulations presented in this 
paper. 



2.1 Assumptions 

The simulations have been performed assuming a time step 
of At = 0.1 Gyr. This quantity has been chosen to be larger 
than the cluster crossing time, in order to ensure that all the 
stars of a given mass group can approach the radial distribu- 
tion predicted by their own King profile, and short enough 
to avoid significant variations of the cluster structural pa- 
rameters. During each time step, the equivalent effect of the 
multiple interactions that occur for a binary population is 
calculated. 

The width of the layers has been chosen to allow a 
higher resolution in the inner (i.e. more populated) regions 
of the cluster. I assumed Ar = 0.1 rc for r < 5 r^, Ar = rc 
for 5 < r/rc < 50 and Ar = 10 rc for r > 50 rc. 

The cluster stars have been divided in mass bins of size 
0.1 Mq in the mass range 0.1 < M/Mq < 5. An additional 
bin formed by stars in the mass range 5 < M/Mq < 120 has 
been considered. These massive stars have lifetimes shorter 
than the time-step of the simulation and are removed during 
the first stage of evolution (see Sect. [2]). I adopted the Ini- 
tial Mass Function (IMF) by Kroupa (2002), which can be 



written as a broken power law dN — m" dm, with a = —1.3 
for 0.1 < M/Mq < 0.5, and a = -2.3 for M > 0.5Mq. 

The initial mass-ratios and periods distributions of the 
binary population have been chosen in agreement with the 
results of Duquennoy & Mayor (1991). In particular, the 
initial mass-ratios distribution has been derived by random 
associating stars belonging to different mass bins. I assumed 
a log-normal distribution of periods g{log P) centered at 
log P{d) = 4.8 and with a dispersion aiog p — 2.3, where 
P is the period expressed in days. The correspondig semi- 
axes distribution of each group of binaries can be derived by 
considering the third Kepler's law 



G (mi -fma) P" 
47r2 



1 1/3 



The distribution of periods has been truncated assuming 
an upper limi10 of P = lO'^ d in agreement with Sills et 
al. (2003). Moreover, a lower limit has been imposed by re- 
jecting all the binaries whose major semi-axes (a) turn out 
to be lower than their limiting orbital distance (fflmin,; see 
Sect. 12.6 (1. The distribution of binary eccentricities has been 
assumed to follow a thermal distribution with probability 
density p{e) — 2e. 

The mass-radius relation and the evolutionary 
timescales of each mass group have been derived from the 
models by Pietrinferni et al. (2006) assuming a metal mass 
fraction Z = 10"^. 

In order to limit the number of free parameters, I as- 
sumed the King parameter Wo as a linear combination of 
the central density and velocity dispersion 



Wo 



log 



Po 



M, 



pc 



1.921+2.934(2) 



This relation has been derived empirically using the cen- 
tral densities and velocity dispersions reported by Djorgov- 
ski (1993) and the cluster concentrations by Trager, King & 
Djorgovski (1995). In the same way, the core radius rc has 
been assumed to be linked to the above parameters in the 
following way 



\Km 



Mc 



102-67Mp 



-1/2 



pc Bellazzini et al. (1998) 



The dependence of the obtained results on the choice of 
some of these assumptions is investigated in Sect. [31 



2.2 Tidal capture 

A mechanism proposed to form binary systems is based on 
the tidal capture of a companion induced by a close en- 
counter between two single stars (Clark 1975). Indeed, when 
two colliding stars approach with a relative velocity and an 
impact parameter lower than given limits, they can form a 
bound system. The rate of binary systems formed through 
this mechanism has been evaluated using the cross section 
(Ttc provided by Kim & Lee (1999). I considered the colli- 
sions between two non-degenerate stars with a polytropic 



^ I did not assumed a criterion to link the upper truncation of 
the period distribution to the local velocity dispersion, assuming 
that binaries born in all clusters with the same initial properties. 
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index x = 1.5 (see their eq. 6). In a time interval At, the 
number ol tidal captured binaries having primary and sec- 
ondary masses mi and m2 and period P turns out to be 



ANb{mi,m2,log P) = 



dNt^tc ,3 



where 



dNt.tc _ ns,r{mi) Ws,r(m2) At gjlog P) 
dV ^ l + 5i,2 



(1 - FcolLv)v (Ttc f{v, r) d V 



where 5i,2 is the Kronecker delta which is unity if stellar 
type 1 and 2 are identical and zero otherwise. The quantity 
Fcoii,v is the fraction of collisions which lead to the formation 
of an unstable binary system whose natural evolution is the 
coalescence in a single massive star (see Sect. 12. 6p and can 
be written as 



FcoU 



As can be seen from the above equations, the binary systems 
formed via tidal capture are assumed to follow the original 
period distribution. 



2.3 Collisional hardening and binaries ionization 

When a single star collides with a binary systems, the two 
objects exchange part of their energy. In particular, the col- 
liding star gain (lose) a fraction of its kinetic energy, which 
is balanced by an decrease (increase) of the binary system 
binding energy, if the hardness parameter 



where 



Et = 



G mi m2 
2a 



(3) 



is larger (smaller) than unity (Heggie 1975). The vari- 
ation of binding energy of the binary system translates into 
a change in the period (i.e. in the orbital separation). This 
process is called "collisional hardening (softening)". If the 
binding energy of the binary system drop to zero then the 
binary is disrupted (binary ionization). 

The code follows the evolution of the binding energy of 
each sub-population of binaries and calculates the number of 
binaries which undergo ionization. In particular, the fraction 
of binding energy gained (lost) by the binary system after a 
time interval At has been assumed to be 



A / , r.^ / i dEb\ 

At J nt.r{mi,m2,logP) y——^ j 



and 



/ 1 dEb\ 
\Eb dt ), 



7.6^/7T G a <ns > if 77 > 1 
'^v .r 

16\/7r G a <nq > loqA .r -1 

■< — n — ^-j — if T7 < 1 



(Binney & Tremaine 1987; Hills 1992). 

Here the quantities < X > indicates the mean value of X 
measured at the radius r. I assumed A = I/77, which is 
a suitable choice for r) < 1 (Binney & Tremaine 1987). I 
neglected the change of the binary eccentricity distribution 
due to coUisions. 



Collisions between binaries plays also an important role, 
in particular when high primordial binary fractions are con- 
sidered. Binary-binary interactions are very efficient for ion- 
izing one (or both) of the two binaries (Mikkola 1983a; 
Sweatman 2007). In a collision between binaries the most 
probable outcome is the disruption of the softer binary and 
hardening of the harder one. As a first-order approxima- 
tion, the collisions between two binary systems have been 
treated as normal collisions between a binary and a single 
point-mass object with a mass equal to the systemic mass 
of the colliding binary system (mi -I- m2). Note that this 
simplification underestimates the amount of energy trans- 
ferred between the two binaries since the geometrical part 
of the cross sections defined above would be few times larger 
then what assumed. Moreover, under this assumption only 
one of the various possible outcomes of a binary-binary colli- 
sion is considered (see Fregeau et al. 2004). However, despite 
the adopted simplifications, the treatment described above 
accounts for most of the contribution of collisions between 
binaries to the ionization process. 

Similarly, the variation of periods and semi-axes have 
been calculated considering that 

^ = and AlogP=-^-^^ (4) 

a Eb 3 Eb 

2.4 Exchanges 

A possible occurrence in a collision between a single star 
and a binary system consists in the formation of an unsta- 
ble triple system which naturally evolves by expelling one 
of the components (usually the least massive) of the origi- 
nal system. The number of exchanges has been estimated by 
means of the cross section Oexch provided by Heggie, Hut & 
McMillan (1996) (see their eq. 4 and 17). Since this formula- 
tion is valid only in the case 7; > 1, I calculated the number 
of exchanges only in this regime. On the other hand, the pro- 
cess of exchange in soft binaries is largely less efficient than 
the one of ionization and can be neglected. Therefore, the 
number of exchanges involving an hard binary having com- 
ponents of masses mi and m2 and period P with a colliding 
star with mass m. is 



dNexch ,3 

; a, r 

dr 



with 



dNexc 

dr 



r{ms) nb^r{mi,m2,log P) Atx 



V cTexch fiv,r) d'^v 



At each time-step, the number of binary sistems and single 
stars have been corrected for this effect. 



2.5 Stellar evolution 

As time passes, most stars evolve and die after timescales 
which depend mainly on the star's mass. Massive stars (M > 
7 Mq) at the end of their evolution explodes as SNII. Less 
massive stars suffers strong mass losses during their latest 
stages of evolution and terminate their evolution as white 
dwarf. 
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When a star belonging to a binary system evolves there 
are many possible outcomes: If the mass of the primary star 
is large enough to produce a SNII explosion part of the 
elapsed energy is converted in a velocity kick that increase 
the total energy of the system. If the total energy (potential 
and kinetic) is positive then the system is disrupted and its 
components will evolve separately. Else, the system remains 
in a bound state in which one of the two companions is the 
compact remnant of the evolved star. 

Otherwise, if the primary component of the binary sys- 
tem is small enough to avoid the ignition of the triple-a 
cycle, mass losses produce a rapid halting of the nuclear 
reactions leading to the formation of a white dwarf (WD). 
Given the large number of low-mass stars, this event is the 
most probable among the various possible outcomes. 

In particular, I considered two important cases: i) the 
envelope of the evolving stars reach the Roche-lobe of the 
companion, igniting the process of mass-transfer (see Sect. 
I2.6|l : a) the star follows an unperturbed evolution leading 
to the formation of a binary system formed by a WD and a 
companion. 

To account for these processes, at each time-step of the 
simulation, a limiting mass has been associated. Then, For 
each mass bin, the code calculates the fraction of stars with 
masses larger then this limiting mass, assuming the stars 
to populate uniformly the mass bin. When this process take 
place in a binary system, the occurrence of the various possi- 
ble processes described above has been considered. In partic- 
ular, in the case of a massive primary component I followed 
the prescriptions of Belczynski et al. (2002). In the case of 
binary systems formed by low-mass stars a distinction be- 
tween processes i ) and ii ) has been made using the criterion 
described in Sect. 12.61 The amount of mass lost at the end 
of the evolution has been calculated using the prescriptions 
by Weidemann (2000). 

The number of binary and single stars is updated ac- 
cordingly. 

2.6 Direct collisions and Coalescence 

From the early '50 it is well known that in globular clusters 
exists a population of massive objects that, in the color- 
magnitude diagram, lie along an extension of the Main Se- 
quence, in a region which is brighter and bluer than the turn- 
off (Piotto et al. 2004 and references therein). These stars 
are called Blue Straggler Stars (ESS). Two mechanisms have 
been proposed to explain the origin of these stars: I) direct 
collision between two single stars (Hills & Day 1976) and 
II) mass-transfer activity in close binary systems (McCrea 
1964). 

Case I) occurrs if the impact parameter of the collision 
between two single stars is small enough to form an unstable 
system in which the process of mass-transfer leads to the 
formation of a single massive star. 

In case II) when the orbital separation of a binary sys- 
tem reduces to a critical distance, the binary system became 
unstable and the ignition of the process of stable mass- 
transfer occurs. This phenomenon can be induced by two 
different mechanisms: Ila) coUisional hardening (see Sect. 
I2.3|l and lib) off-Main Sequence evolution of the primary 
component of the binary system. In the former case, hard 
binaries subject to a high number of collisions tend to re- 



duce their orbital separation (see eqn.|4]). In the latter case, 
when the primary component of the binary system exhausts 
the hydrogen in its core, it expands its envelope by a large 
factor (which depends on its mass) possibly reaching the 
Roche-radius of the secondary component. 

The critical distance between two companions in a bi- 
nary system to ignite stable mass-transfer is 



= 2.17 



mi + 1712 
m2 



1/3 



7?2 Lee & Nelson (1988) 



The final product of these processes is a system contain- 
ing a massive star having a mass mi < m < (mi +m2). The 
formation of BSSs afi'ects the number of binary and single 
stars. 

To evaluate the frequency of the processes described 
above, I first consider the cross-section acoii for a stellar 
collision which leads to a close binary system with orbital 
separation a < amin (case 7). This cross-section has been 
calculated using the formulation by Lee & Ostriker (1986; 
see their eqns. 2.8 and 4.12) 

2Gm2 

^coll — T^dmin ^ 

Then, the number of coalesced stars via direct collision 
can be written as 



ANt{mi,m2,log P) = 



dNi,^coii{mi,m2) ,3 
— a r 



where 

dNi^coH _ ns,r(mi) ns,r{m2) At 



dr 



l + 5l,: 



dr 



V (Tcoii f{v,r) d V 



The fraction of coalesced stars via coUisional hardening 
(case Ila) has been calculated by following the evolution of 
the semi-axes of the binary population (see Sect. 12. 3|) . In- 
stead, when the coalescence is due to the evolution of the 
primary star of a binary system (case lib), the code calcu- 
lates the fraction of system with an evolving primary star 
which satisfies the condition 



a < a„ 



R\,'max 

^ Ri 



Where R\ and Ri,max are the radii of the primary star dur- 
ing its quiescient and evoluted stages, respectively. 



2.7 Evaporation 

During its evolution, a globular cluster undertakes the tidal 
stress of its host galaxy. As a result of this interaction, it 
lose part of its stellar systems decreasing its total mass. The 
rate of evaporating systems depends on the cluster orbit and 
on its structural parameters (Gnedin & Ostriker 1999). 

Here, the treatment of evaporation has been greatly 
simplified assuming the cluster to lose a constant fraction 
V of its objects during all its evolution. Although this is a 
crude approximation, it can be used to simulate the effect 
of a tidal field on the relative amount of single and binariy 
stars survived to the evaporation process. A homogeneous 
evaporation rate is also in agreement with the most recent 
N-body simulation (Hurley et al. 2007; Kim et al. 2008). At 
a given distance from the cluster center, I assumed the ve- 
locity distribution of a population of objects with mass m 
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can be well represented by a Maxwellian distributiorlfl 

where K is a term proportional to and A is a normal- 
ization factor. The quantity K depends on the stellar mass 
according to the status of relaxation of the cluster (see Ap- 
pendix A) . 

The relative number of evaporating systems belong- 
ing to a given group of mass can be therefore calculated 
by integrating over the entire cluster extension the frac- 
tion of stars with a velocity higher than the escape velocity 



Jo 



dr 



where 



du' (m) 
dr 



1/2 -u I 

y e 0-y 



l3Wr(m) 



Wr{m) = 



The term P accounts for the state of partial relaxation (see 
Appendix A) 



where 



[7 < m > -|-(1 — ■y)m] 

Nm{r < Vrel) 



7 ^ 



is the fraction of stars with mass m located in the region 
where relaxation already occurred. 

Once this ratios are calculated for all the mass groups, 
the normalization factor B can be derived by imposing that 



/\Nsys = ^ iV'(mi) = vNs. 



2.8 Cluster dynamical evolution 

During its evolution, the cluster core lose part of its mass as 
a result of the processes of mass-segregation and evapora- 
tion. To mantain its virial equilibrium the cluster collapses 
increasing its central density and velocity dispersion. A sig- 
nificant population of hard binaries feeds kinetic energy into 
the cluster through binary-single and binary-binary interac- 
tions which counteract the effects of evaporation. This pro- 
cess can halt and even reverse the contraction of the core 
(Hills 1975; Gao et al. 1991). Fregeau et al. (2003) showed 
that also a small fraction of hard binaries is sufficient to 
support the core against collapse significantly beyond the 
normal core-collapse time. 

I considered the evolution of the structural properties of 
the cluster in a highly simplified way. In particular, the code 
calculates at each time-step of the simulation the number 



of stars contained in the cluster core and updates the core 
radius and the central velocity dispersion using the relations 

and 



'eJT' 



where 



EL 

E 



ln(N'/2) 
ln{N/2) 



AH) 



Hills (1975) 

Where the symbols X and X' indicate the values of X in two 
subsequent time-steps. The term E' /E represents the ratio 
between the cluster binding energies in two subsequent time- 
steps due to the contribution of binaries. Here, all quantities 
refers to the cluster core and ^ refers to the fraction of hard 
binaries. The exponent z((^) is related to the heating of bina- 
ries due to binary-single stars and binary-binary interaction 
and can be written as 

z{0 = (^/l-3) + 10.3 ^-{1.2 - 1/ < 77 >) Mikkola (1983b) 

Moreover, the King parameter has been updated using eq. 
[21 In this way, the cluster structural parameters satisfy the 
virial theorem during all the simulation. 

In Fig. [2] the evolution of the core radius for the mod- 
els s5p3e5f50 and s5p3e5fl0 are shown. As can be seen, the 
general trend of the core radus evolution is characterized by 
a progressive contraction of the core induced by the ever 
continuing losses of stars. In model s5p3e5f50 the heat gen- 
erated by hard binaries halt the contraction and lead to a 
stabilization toward an equilibrium value. The behaviour of 
the evolution of the cluster characteristic radii qualitatively 
agree quite well with the most recent N-body simulation 
computations (see Fregeau et al. 2003; Kim et al. 2008). 
The details of the core oscillations and post-collapse evolu- 
tion are not reproduced by the code as a consequence of the 
simplified treatment described above. However, this should 
have only a minor effect on the predicted evolution of the 
properties of the binary population. 



3 RESULTS AND DEPENDENCE ON THE 
ASSUMPTIONS 

I run a set of 36 simulations spanning a wide range in density, 
velocity dispersion, evaporation efficiency and initial binary 
fraction. A list of all the performed simulations and their 
initial and final parameters is provided in Table 1. 

To test the dependence of the obtained results on the 
adopted assumptions in the following sub-sections I compare 
the model s5p3e5f50 with a set of simulations preformed 
with the same initial conditions but different assumptions. 
The results of these additional simulations are listed in Table 
2. 



^ Although this assumption holds only in the core, I extended its 
validity to the entire cluster only to calculate the relative fraction 
of evaporating stars belonging to different mass groups. 



3.1 Dependence on Metallicity 

As described in Sect. 12.11 the simulations have been per- 
formed adopting a metal mass fraction Z — 10~^. Metal- 
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Figure 1. The evolution of the binary fraction is shown as a function of time for the reference model s5p3e5f50 (solid lines) and for 
the models calculated under different assumptions. Left panels show the binary evolution in the whole cluster, right panels show the 
binary evolution in the cluster core. 



licity influences both the evolutionary timescales and the 
mass-radius relation and can in principle affect the evolu- 
tion of the binary population. 

In the upper panels of Fig. [T] the behaviour of the bi- 
nary fraction as a function of time is showed for the model 
s5p3e5f50 and for a simulation performed with Z — 10~*. 
As can be noted, the two models are indistinguishable even 
when the core binary fraction is considered. Also the final 
binary fraction of the metal-poor model is the consistent 
with the reference metal-rich one within 0.3% (see Table 2). 

I conclude that the metallicity plays a neglegible role in 
the evolution of the binary fraction in the presented simula- 
tions. Note however that metallicity can play an important 
role in the primordial fraction of binaries and in the for- 
mation of X-ray binaries (see Ivanova 2006) . In the models 
presented here such effects can remain hidden because of the 
adopted simplifications. 



3.2 Dependence on the IMF 

Another important ingredient of the simulations is the IMF. 
Indeed, it determines the relative fraction of each mass group 
and influences the distribution of mass ratios. 

Although the power-law exponent of the IMF is well 
estabilished for stars with masses M > 0.5Mq, there are 



still large uncertainties in the shape of the low-mass end 
of the IMF. Moreover, some authors claim a dependence 
of the power-law exponent of the low-mass end of the IMF 
on metallicity (McClure et al. 1986; Djorgovski, Piotto & 
Capaccioli 1993). 

In the central panels of Fig. [T] the model s5p3e5f50 is 
compared with two simulations performed using two ex- 
treme values of the IMF power-law exponent in the mass 
range M < 0.5Mq {a = and a = 1.8). As can be seen, 
the three models show a very similar behaviour during all 
the cluster evolution. The flnal binary fractions of the three 
models are also compatible within ~ 5% . 

Therefore, although the IMF plays a fundamental role 
in the determinations of the main characteristics of the bi- 
nary population, it plays only a minor role in the evolution 
of the binary fraction. 



3.3 Dependence on the Period distribution 

One of the most important assumptions made in the sim- 
ulations regards the distribution of periods of the binary 
population. The shape and the extremes of this distribution 
determine the ratio of soft-to-hard binaries (i.e. the fraction 
of binaries which can undergo ionization) . 

To evaluate the impact of these assumptions on the evo- 
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Table 1. Complete list of simulations. 



initial conditions final conditions 
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10 





8.3 
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10.0 
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Table 2. Results of the test simulations. 



initial conditions final conditions 



Model ID 
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% 
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s5p3e5f50 


50 


-3 


1.3 


gaussian 
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29.9 


39.0 


s5p3e5f50.zlow 


50 


-4 


1.3 


gaussian 


7 


30.0 


39.3 


s5p3e5f50.ml8 


50 


-3 


1.8 


gaussian 


7 


24.9 


35.7 


s5p3e5f50.m0 


50 


-3 


0.0 


gaussian 


7 


31.3 


41.2 


s5p3e5f50.pflat 


50 


-3 


1.3 


fiat 


7 


39.8 


48.3 


s5p3e5f50.p6 


50 


-3 


1.3 


gaussian 


6 


35.2 


44.9 



lution of the binary fraction I compared the model s5p3e5f50 
with two simulations performed assuming i) a flat period 
distribution truncated at log P{d) = 7 (where P is ex- 
pressed in days) and ii) a. log- normal distribution truncated 
at log P{d) — 6. The results of this comparison are shown 
in the lower panels of Fig. [T] Note that, while a variation 
of the upper limit of the period distribution produces only 



small changes (A^ < 6%) in the final binary fraction, the 
model with a flat period distribution mantains a significantly 
higher fraction of binary systems during all the cluster evo- 
lution. In particular, after 13 Gyr, this last model contains 
^ 10% of binaries more than the reference model. As pointed 
out above, this effect is due to the higher fraction of hard 
binaries resulting from the flat distribution of periods (60% 
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Figure 2. The evolution of the core radius is shown as a func- 
tion of time for the models s5p3e5f50 (solid line) and s5p3e5fl0 
(dashed line). 




Figure 3. Derivative of the binary fraction as a function of time 
for model s5p3c5f50. The contributions of the various processes 
are indicated. 



at the beginning of the simulation) with respect to the log- 
normal distribution (43%). 

All the results presented in the following sections are 
based on simulations run assuming a log-normal distribu- 
tion of periods truncated at log P{d) = 7. Although the 
dependence of the properties of the binary population on 
the cluster structural and environmental parameters should 
not be influenced by this choice, it is important to bear in 
mind that the absolute results of the individual simulations 
are highly sensitive to this assumption. 



4 EVOLUTION OF THE BINARY FRACTION 

The simulations listed in Table 1 cover a range of structural 
and environmental parameters comparable to those spanned 
by Galactic globular clusters. Thus, they can be used to 
study the efficiencies of the various processes of binary for- 
mation/destruction as a function of the initial conditions. 
For this purpose consider the derivative 

dtg ^ \dtgJx 

This quantity can be seen as the sum of the contributions 
of the individual processes, where 

/ d£\ _ NsysdxNb - NtdxNsys 
\dtgJx N'iys 

represents the contribution of the process X to the above 
derivative. In Fig. |3] the behaviour of the derivatives de- 
fined above are shown as a function of time for the model 
s5p3e5f50. Such a behaviour is qualitatively the same in 
the other models with different initial conditions. As can be 
noted, the process of ionization is the main responsible for 
the variation of the binary fraction reaching its maximum 



efficiency in the first Gyrs of evolution. Another important 
contribution is given by evaporation (in models where i/ ^ 0) 
which enhances the fraction of binaries sistematically dur- 
ing the whole cluster evolution. All the other processes play 
only a minor role. 

In Fig. 13] and [S] the behaviour of the binary fraction as 
a function of time is shown for all the simulations with ini- 
tial binary fraction ^ — 50%. As can be noted, the general 
evolution of the binary fraction is characterized by a sud- 
den decrease in the first Gyr, followed by a constant trend 
during the subsequent Gyrs of evolution. In simulations run 
assuming i/ 0, the fraction of binaries slightly increase in 
the last Gyrs of evolution. The fraction of binaries in the 
core is sistematically larger than that of the entire cluster. 
This difference is more evident after 2-3 Gyr evolution. 

The strong initial decrease of the binary fraction is due 
to the process of binary ionization which destroy a large 
fraction of binaries during the first Gyr of evolution. After 
this initial stage, all the soft binaries have been destroyed 
and the evolution of the binary fraction is driven by the other 
processes which have a smaller efficiency. In particular, if the 
efficiency of the process of evaporation is large, the fraction 
of binaries increases again. Binary systems are in fact on 
average more massive then single stars and tends to have 
therefore smaller velocities. Therefore, a larger number of 
single stars reach the cluster escape velocity with respect 
to binary systems, thus increasing the relative fraction of 
binaries. As relaxation proceeds a larger fraction of binaries 
sink into the central regions of the cluster, as a consequence 
of the process of mass segregation, thus increasing the core 
binary fraction. 

Fig. |6] shows a map of the fraction of surviving binary 
systems after 13 Gyr in the log po — <yv,o plane for the sim- 
ulations with initial binary fraction ^ — 50%. As can be 
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Figure 4. The evolution of the binary fraction in the core (dashed lines) and in the whole cluster (solid lines) are shown as a function 
of time for all the models with 1/ = and f = 50%. 



seen, a clear trend is visible as a function of the cluster final 
conditions. In particular: 

• The fraction of binaries decreases by increasing the cen- 
tral velocity dispersion; 

• The fraction of binaries decreases by increasing the cen- 
tral density; 

• The fraction of binaries increases by increasing the ef- 
ficiency of the evaporation. 

The physical reasons at the basis of these trends are 
linked to the efficiencies of the processes of binary ionization 
and evaporation. Indeed, by increasing the cluster velocity 
dispersion increase the ratio of soft-to-hard binaries (see eq. 
[3| which are destroyed by the process of ionization in the 
first Gyrs of evolution. The efficiency of ionization increase 
also by increasing the cluster density (see eg. I2.3|l . On the 
other hand, evaporation tends to increase the fraction of 
binaries (see above). 

Interestingly, the efficency of ionization and the clus- 
ter mass have the same dependences on the cluster struc- 
tural parameters (po and a^fi)- Indeed, more massive clus- 
ters have, on average, larger densities and velocity disper- 
sions (Djorgovski & Meylan 1993). In Fig. [7] the fraction of 
binaries after 13 Gyr is plotted against the residual cluster 
mass for all models with ^ — 50% and u — 5%. As expected, 
a clear anticorrelation is clearly visible between these two 



quantities. The same result can be achieved considering dif- 
ferent values of the primordial binary fraction and evapora- 
tion efficency. 

In Fig. [S] the evolution of the binary fraction for the 
model s5p3e5f50 is compared with that calculated for a 
model with the same initial conditions but a significatly 
smaller initial binary fraction (model s5p3e5fl0) . The binary 
fractions for the two models have been normalized to their 
initial values in order to compare the relative trend of their 
evolution. Note that the model which starts with a lower bi- 
nary fraction mantain a larger fraction of binaries in the first 
stages of its evolution. This evidence is more evident when 
the core binary fraction is considered. This is a consequence 
of two effects: i) the effect of the binary-binary collisions and 
ii) the dependence of the local relaxation time on the mean 
stellar mass (see eq. [l|. Indeed, in clusters with higher bi- 
nary fractions, binary-binary collisions are ~ 25 times more 
frequent and, given their higher energetic budget, dominate 
the process of binary ionization. Moreover, the mean stellar 
mass is smaller in clusters with smaller binary fractions. The 
local relaxation time in these clusters is therefore longer, and 
a smaller fraction of binaries is therefore involved in the pro- 
cess of binaries destruction. This effect is amplified in the 
core where more binaries sink as a result of mass segrega- 
tion. The trend is reversed after few Gyrs when the impact 
of evaporation became more important. In fact, models with 
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Figure 7. The surviving fraction of binaries after 13 Gyr is shown 
as a function of the residual cluster mass for all the models with 
^ = 50% and u = 5%. A clear anticorrelation between these 
quantities is visible. 

higher fractions of binaries tend to mantain binaries more 
efRcently. At the end of the simulation the two models con- 
tain a similar fraction of their initial budget of binaries in 
the core regardless of their primordial content. In the follow- 
ing I discuss the dependence of the evolution of the binary 
population considering the models calculated assuming an 
initial binary fraction ^ = 50%. The conclusions are qualita- 
tively the same for the models with a smaller initial binary 
fraction. 



5 RADIAL DISTRIBUTION OF BINARIES 

Binary systems are on average more massive then single 
stars. Therefore, they are expected to have a different ra- 
dial distribution with respect to single stars, as a result of 
the process of mass segregation. 

Fig. [5] shows the fraction of binaries (model s5p3e5f50) 
normalized to the central value as a function of the radial 
distance from the cluster center at different times. As ex- 
pected, the fraction of binaries shows a central peak, de- 
creasing toward the external regions of the cluster. As time 
passes, relaxation occurs in the external regions and the 
cluster increases its concentration as a result of the ever 
continuing losses of stars. Consequently, the fraction of bi- 
naries populating the outermost part of the cluster slightly 
increase. 

In Fig. [To] the radial behaviour of the normalized frac- 
tion of binaries (model s5p3e5f50) calculated after 13 Gyr 
is compared with models with extremely different initial 
conditions (models s2ple5f50 and sl3p5e5f50; upper panel) 
and with a smaller evaporation efficiency (s5p3e0f50; bottom 
panel). Note that the slope of the decreasing trend of the bi- 
nary fraction in the external part of the cluster is less steep 
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Figure 8. The binary fraction normalized to the initial value is 
shown as a function of time for the model s5p3e5f50 (solid lines) 
and s5p3e5fl0 (dashed lines). The upper panel refers to the entire 
cluster, the lower panel to the cluster core. 

in model sl3p5e5f50 which undergo a stronger dynamical 
evolution. The same behaviour is visible when evaporation 
accelerates the cluster dynamical evolution. 

Summarizing, the slope of the external decrease of the 
radial distribution of the binary fraction is steeper as the 
cluster is dynamically younger. 

A note of caution is worth: the code assumes that all 
the cluster stars follow the radial distribution predicted by 
a given King model according to their masses during the 
whole cluster evolution. Under this assumption, the process 
of mass segregation removes instantaneously all the radial 
differences in binary fraction produced by the various mech- 
anisms of binary formation/destruction. However, while this 
approximation is reasonable in the central regions (where 
the local relaxation time is relatively short) in the outer 
regions this assumption could not hold. Therefore, the pre- 
dicted radial distribution of binaries could not be reliable in 
the external regions of the cluster. 



6 PERIOD DISTRIBUTION 

As outlined in Sect.|3J one of the main processes that drives 
the evolution of the binary fraction is the process of binary 
ionization. The binaries which are more subject to this pro- 
cess are those with smaller binding energy (i.e. longer pe- 
riods). The shape of the distribution of periods therefore 
changes during the cluster evolution according to the effi- 
ciency of the process of binaries ionization. 

In Fig. [TT] the initial and final periods distributions of 
three models with different ionization efficencies are shown. 
As expected, in models with a high ionization efficency (see 
model sl3p5e5f50) all binaries with log P{d) > 2 are de- 
stroyed. On the oppsite side, when ionization plays a minor 
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Figure 5. Same of Fig.Hbut for all the models with u = 5%Gyr-'^ and g = 50%. 



role (model s2ple5f50) a significant fraction of binaries with 
period as long as log P{d,) = 6 still survives. 



7 MASS-RATIO DISTRIBUTION 

During the cluster evolution also the mass-ratios distribu- 
tion changes. There are three main processes that drive its 
evolution: stellar evolution, ionization and exchanges. In- 
deed, at the beginning of the simulation a large number of 
binaries are formed by a massive primary star which dies 
during the cluster evolution. The minimum possible mass- 
ratio therefore increase with time. Moreover, the mean bind- 
ing energy of the binary population is proportional to the 
average mass-ratio. Thus, in clusters where ionization has a 
higher efficiency the mass-ratios distribution should be more 
efficiently depleted up to larger values of q. During close en- 
counters between a binary system and a massive colliding 
single star, the secondary star of the binary system is usu- 
ally ejected from the system. When this process becomes fre- 
quent, the more stable systems against exchanges are those 
formed by equal mass components. Thus, clusters with high 
coUisional rates should show an increase of the number of 
equal mass binaries. 

In Fig. [T5] the distribution mass-ratios of the model 
s5p3f5e50 is compared with those of two models with ex- 
tremly different densities (models s5p2e5f50 and s5p5e5f50) 



and evaporation efficiency (model s5p3e0f50). As expected, 
in all the distributions there are few systems with q < 0.2, as 
a result of stellar evolution. Moreover, the model s5p5e5f50 
shows a lack of stars with low mass-ratios with respect to 
the three other models. This is an effect of the high efficiency 
of ionization which depletes the distribution of mass ratios 
up to larger values of q. In none of the performed simula- 
tions an increase of the fraction of equal mass binaries is 
noticeable. This means that the efficiency of the exchange 
process is rather low over the entire range of parameters 
spanned by the simulations. Indeed, even in the densest sys- 
tems, where collisions are more frequent, the fraction of ex- 
changes Nexch/Nb never exceeds few percent. 



8 COMPARISON WITH N-BODY AND 
MONTE CARLO SIMULATIONS 

The analytical model presented here has been compared 
with the most recent existing N-body (Hurley et al. 2007) 
and Monte Carlo (Ivanova et al. 2005) simulations. As out- 
lined in Sect. [T] the two above approaches are completely 
different and lead to apparently discrepant results. In fact, 
while Monte Carlo simulations by Ivanova et al. (2005) pre- 
dict a strong depletion of the core binary fraction, N-body 
simulations by Hurley et al. (2007) show the inverse trend. 
Part of this discrepance can be due to the very different 
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Figure 6. Maps of the surviving binary fraction in the plane log po — crv,o for ^ = 50%, u = {left panels) and u = 5%Gyr~^ 
{right panels). Lower panels refers to the core binary fraction, upper panels to the whole cluster binary fraction. Darker regions 
indicate higher binary fractions in steps of 10%. The locations of the Galactic globular clusters (from Djorgovski 1993) are also marked 
with open circles. 



initial conditions of the two simulations (Fregeau 2007). In- 
deed, Ivanova et al. (2005) simulated a series of high density 
cluster with a significant velocity dispersion. In the follow- 
ing I refer to their model BOS which has n = 10^ pc'^ and 
o'v,o = 10 Km . In the approach followed by these au- 
thors, the evolution of the binary population is supposed to 
proceed in a "fixed backgroung" (i.e. a cluster whose struc- 
tural parameters do not change during the simulation). The 
fraction of evaporating stars in these simulation is also very 
small (~ 6% after 13 Gyr; N. Ivanova, private communi- 
cation). This condition can be reproduced by assuming a 
very small efiicency of evaporation (i^ ~ 0). Instead, Hurley 
et al. (2007) performed four different simulations assuming 
lower densities and velocity dispersions for cluster which in- 
teract with the Galactic disk on a circular orbit in the solar 
neighborhood. I considered their model K24-50 which have 
an initial binary fraction ^ = 50% a central star density 
n = 100 pc~^ and a velocity dispersion of (t„,o = 3 Km . 
According to Fig. 4 of Hurley et al. (2007), the cluster lose 
~ 70% of its object after 13 Gyr (comparable with a model 
with V = 8%Gyr~^). Considering the results obtained in 
Sect. |4] (see e.g. Fig.[6l the Hurley et al.'s model is expected 
to retain much more binaries with respect to that of Ivanova 
et al. (2005). 



Moreover, the two above authors make different as- 
sumptions regarding the distribution of periods. In fact, al- 
though both authors adopt a flat distribution, Ivanova et 
al. (2005) assume an upper truncation at log P{d) = 7 
while Hurley et al. (2007) impose an upper limit to the or- 
bital distance of 50 AU (corresponding to a truncation at 
log P{d) ~ 5, see their Fig. 2). Consequently, the Ivanova et 
al.'s model contains a larger initial fraction of soft binaries 
which are easily destroyed during the cluster evolution. 

To test the consistency of the approach followed in this 
paper with the ones quoted above, I run two simulations 
with the same assumptions and initial conditions of the 
above authors. I will refer to these two simulations as the 
H-like and I-like models. Considering that both simulations 
start with an initial binary fraction ^ = 50%, at the end of 
its evolution the H-like model contains a fraction of binaries 
in the core of ^ = 76.9% while the I-like contains only a 
fraction of binaries f = 25.1%. Although the absolute val- 
ues of the above simulations differs from those reported by 
those authors, the qualitative behaviour of both simulations 
is well reproduced. This means that all the three approaches 
are qualitatively self-consistent, and that the observed dis- 
crepancy (A^ ~ 52%^ is due to both the intrinsic different 
starting parameters and the different assumptions. 
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Figure 9. Normalized binary fraction as a function of the radial 
distance to the cluster center at different times from the beginning 
of the simulation for the model s5p3e5f50. 



r/r 

5 ' 10 




0.4 - 



o.a - 



5 10 

Figure 10. In the upper panel the comparison among the nor- 
malized binary fractions as a function of the radial distance 
to the cluster center is shown for the model s5p3e5f50 (solid 
line), s2ple5f50 (dashed line) and sl3p5e5f50 (dotted line). The 
lower panel shows the same comparison but for the models 
s5p3e5f50 (solid line) and s5p3c0f50 (solid line). 



To quantify the effect of the different assumptions on 
the resulting binary fraction, I performed a series of sim- 
ulations adopting the starting conditions of the H-like and 
I-hke models but different values of the maximum period 
and evaporation efficiency. The complete set of simulations 
are summarized in Table 3. As can be seen, adopting the 
same assumptions, the difference between the fractions of 
survived binaries predicted by the two models reduces by 
a factor of two (A^ = 21 27%). Therefore, at least an 
half of the difference between the fractions of survived bi- 
naries predicted by the above models is due to the different 
assumptions. The remaining difference can be addressed to 
the different initial cluster conditions. 



9 COMPARISON WITH OBSERVATIONS 

The code described here has been used to interpret the ob- 
servational evidence that young globular clusters contain a 
higher fraction of binaries with respect to older ones (Sol- 
lima et al. 2007). As outlined in Sect. O the fraction of 
binaries should not show significant variations in the last 
Gyr of evolution. Therefore, this trend could be due to i) 
differences in the primordial binary fractions or ii) to a 
different evolution of the binary populations linked to the 
different structural and environmental parameters of these 
two groups of clusters. Indeed, the group of young globu- 
lar clusters of the sample of Sollima et al. (2007) is char- 
acterized by an age of ~ 7 Gyr, a mean central den- 
sity < logpo >= 0.65 (where po is expressed in Mqpc~^) 
and a velocity dispersion < a,j,o >~ l-3Km which 
are significantly different from those of the old group of 
that sample (< tg >— llGyr, < logpo >= 2.15 and 
< cFvfi >= 3-SKm s~^). According to the results found in 
Sect.m the group of young clusters should mantain a larger 
fraction of binaries during their evolution. 

To test if these differences can be responsible for the 
difference in the observed binary fractions, I tried to repro- 
duce the mean observed conditions of the two above groups 
of globular clusters assuming the same initial binary frac- 
tion. It is worth of noticing that the fraction of binaries 
derived by Sollima et al. (2007) is calculated in a restricted 
range of masses and mass-ratios. To compare the results 
of the simulations with the observations, I compared the 
minimum binary fraction ^min defined in Sollima et al. 
(2007) with the ratio between the number of binaries with 
a primary star in the mass range 0.5 < mi/M© < 0.8 and a 
mass ratio q > 0.5, and the number of objects in the same 
mass range. Of course, the relation between this quantity 
and the global binary fraction depends on the distribution 
of mass-ratios. Given the arbitrary choice of this distribu- 
tion (see Sect. 12. ip the absolute value of the initial frac- 
tion of binaries derived here is not reliable and can be only 
used for differential comparisons. According to Sollima et al. 
(2007), the groups of old and young globular clusters have 
minimum binary fractions (,min = 6% and ^min ~ 20% 
respectively. 

The model that better reproduces the fraction of bi- 
naries in the old group of globular clusters after 11 Gyr 
is characterized by an initial binary fraction of ^ = 8%. I 
found no combinations of the initial parameters able to re- 
produce the observed fraction of binaries in the young group 
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Table 3. Results of the H-like and I-like simulations. 
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Figure 11. Period distributions of three different models. Grey histograms show the distribution after 13 Gyr. Open histograms show 
the distribution at the beginning of the simulation. Periods are expressed in days. 
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Figure 12. The mass-ratios distributions after 13 Gyr of models s5p2e5f50 {left panel) and s5p5e5f50 (right panel) are shown with 
grey histograms. The mass-ratios distribution of model s5p3e5f50 (open histograms) is shown in both panels for comparison. 
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of globular clusters without assuming a larger initial binary 
fraction. This indicates that, although the dynamical status 
of the young group of globular clusters favors the survival of 
a larger fraction of binaries, the most of the observed differ- 
ence have to be due to primordial differences in the cluster 
binary content. 



10 CONCLUSIONS 

I presented a code designed to simulate the evolution of the 
properties of a binary population in a dinamically evolving 
star cluster. A number of simulations spanning a wide range 
of structural and environmental parameters have been run. 

In general, the fraction of binaries appears to decrease 
with time by a factor 1-5, depending on the initial cluster 
parameters. In particular, the fraction of binaries quickly 
decreases in the first Gyrs of evolution, as a result of the 
high efficiency of the ionization process in this initial stage. 

The analysis of the contributions of each mechanism of 
formation and distruction of binaries indicates that the main 
processes that drive the evolution of the binary fraction in 
globular clusters are the processes of binary ionization and 
evaporation. This result seems to be confirmed by the ob- 
servational fact that open clusters and low-density globu- 
lar clusters contain more binary systems than dense high- 
velocity dispersion globular clusters (Sollima et al. 2007). 
As a consequence, the fraction of survived binaries increases 
when the cluster structural parameters support a lower effi- 
ciency of the process of ionization and in clusters subject to 
strong evaporation. Moreover, also the final distributions of 
periods and mass-ratios change as a function of the cluster 
structural parameters. In particular, clusters with a higher 
efficiency of binary ionization tend to have binaries with 
shorter periods and higher mass-ratios. At present, the range 
in central density and velocity dispersion spanned by the 
sample of globular clusters with homogeneous estimates of 
binary fraction is very small and does not allow a proper 
comparison (see Sollima et al. 2007). However, preliminary 
results by Milone et al. (2008) based on the analysis of a 
larger sample of globular clusters seems to confirm this trend 
. These authors found a significant anticorrelation between 
the fraction of binaries and the cluster luminosity (i.e. mass). 
Indeed, more massive clusters have, on average, larger den- 
sities and velocity dispersions (Djorgovski & Meylan 1993) . 
Therefore, the observed anticorrelation could be due by the 
fact that the cluster mass and the efficency of binary ion- 
ization have the same dependence on the cluster structural 
parameters (po and ct^,o). 

The dependence of the obtained results on the assump- 
tions has been tested. The evolution of the binary fraction 
appears to be non-sensitive to the assumed cluster metallic- 
ity and to the shape of the low-mass end of the IMF. A sig- 
nificant dependence on the initial shape of the distribution of 
periods has been found. In particular, a difference of ~ 10% 
has been found by switching from a log-normal to a ffat 
distribution. Therefore, the decreasing trend of the binary 
fraction with time could be even inverted Eissuming a differ- 
ent distribution of periods at least in those clusters where 
the efficiency of the ionization process is small. Unfortu- 
nately, the shape of this distribution is still largely uncertain 
(Halbwachs et al. 2003). Studies performed on large sam- 



ples of binaries in the Galactic field are not able to distin- 
guish between the various proposed distributions, suffering 
strong observational biases (Abt 1983; Eggleton et al. 1989; 
Duquennoy & Major 1991; Halbwachs 2003). Moreover, de- 
tections of long-period binary systems (with P > W^d) arc 
limited by the intrinsic impossibility to detect photomet- 
ric and/or kinematical variability over such large timescales. 
Given the importance of this parameter, an indication of the 
true shape of the periods distribution would be valuable. 

The predicted radial distribution of binary systems 
shows a decreasing trend with central peak and a rapid drop 
toward the external regions of the cluster. The reason at the 
basis of this behaviour is linked to the process of mass seg- 
regation which produces a concentration of the massive bi- 
nary systems in the central region of the cluster. The higher 
concentration of binary systems has been already observed 
by several authors in different clusters (Yan & Reid 1996; 
Rubenstein & Bailyn 1997; Albrow et al. 2001; Bellazzini et 
al. 2002; Zhao & Bailyn 2005; Sollima et al. 2007). 

The results of the code presented here have been com- 
pared with the most recent N-body and Monte Carlo simu- 
lations available in the literature. Despite the many adopted 
simplifications, the predictions of the code appear to be qual- 
itatively consistent with the results of the above approaches. 
I estimated that at least an half of the difference between 
the fractions of survived binaries predicted by the N-body 
simulations by Hurley et al. (2007) and the Monte Carlo 
simulations by Ivanova et al. (2005) is due to the different 
assumptions made by these authors regarding the upper end 
of the periods distribution and the treatment of evaporation. 
The remaining difference can be addressed to the different 
initial cluster conditions assumed by these authors (as al- 
ready suggested by Fregeau 2007). 

The code presented here has been used to interpret the 
differences in the binary fractions measured in the sample of 
globular clusters presented by Sollima et al. (2007). There is 
no combination of initial parameters able to reproduce the 
observed fraction of binaries observed in group of clusters 
with different ages unless assuming a different initial binary 
content. The group of young globular clusters in the sample 
of Sollima et al. (2007) is formed by three clusters (namely 
Terzan 7, Palomar 12 and Arp 2). These clusters are also the 
most distant from the Sun and they are thought to belong 
to the Sagittarius Stream (Bellazzini et. al 2003). Thus, they 
might be stellar systems with intrinsically different origins 
and properties, whose initial conditions could significant dif- 
fer from those of "germine" Galactic globular clusters (see 
also Sollima et al. 2008). Future studies addressed to the 
estimate of the binary fraction in other clusters suspected 
to have an extragalactic origin will help to understand how 
the enviromnental conditions influence the original content 
of binaries of globular clusters. 
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APPENDIX A: PARTIALLY-RELAXED 
MULTI-MASS KING MODELS 

The density and velocity distribution profiles of globular 
clusters are well represented by King models (King 1966). 
These models assume that cluster stars located at a distance 
r from the cluster center follow a distribution of velocities 

f{e,J)=Ae-^'\e-'^ ~1) (Al) 

(Mitchie 1963; King 1966; Gunn & Griffin 1979). 
where K is a quantity proportional to that accounts for 
the cluster relaxation, is the energy for unity of mass 
(er — f&r + i'^/2), J is the angular momentum, and A is 
a normalization factor. The term exp(-/3J^) accounts for 
anisotropies in the distribution function (here I set this 
quantity equal to unity). The shape of distribution of stars 
in these models is complely defined by the King parameter 
Wo = — $o/o"u ■ The density profile can be derived from the 
above distribution considering the Poisson equation and the 
relation 

pr(m) = / f{v,r)cPv 
Jo 

After a timescale comparable to the local relaxation 
time, collisions become frequent and produce the equiparti- 
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tion of the kinetic energy {K oc m~^). Under this condition, 
the distribution of velocities becames a function of the mass. 
Otherwise, the distribution of velocities does not depend on 
mass {K — const) and the above formulation reduces to the 
case of a single-mass King model. However, relaxation oc- 
curs after different timescales at different distances from the 
cluster center (see eq. [1} ■ Therefore, only a fraction of stars 
follow relaxed orbits. 

For this reason, the code calculates at each time-step 
of the simulation the fraction of stars located in the region 
where relaxation already occurred. Defining 

the distribution of velocities has been then assumed to follow 
the form of eg. I All where 

7 < m > +{1 - j)m 2 
m 

The density and velocity dispersion profiles have been there- 
fore calculated accordingly. 



